Stripes due to the next-nearest neighbor exchange in high-T c cuprates 



O 

o 

X> 

PL, 
OS 



Toru Sakai 

Faculty of Science, Himeji Institute of Technology, 3-2-1 Kouto, Kamigori, Ako-gun, Hyogo 678-1297, Japan 

(February 2001) 

We propose a possible mechanism ol the charge stripe order due to the next-nearest neighbor 
exchange interaction J' in the two-dimensional t-J model, based on the concept of the phase sepa- 
ration. We also calculate some hole correlation functions of the finite cluster of the model using the 
numerical diagonalization, to examine the realization of the mechanism. It is also found that the 
next-nearest neighbor hopping t' suppresses the stripe order induced by the present mechanism for 
t' < 0, while it enhances for il > 0. 

PACS Numbers: 71.10.Fd, 71.45.Lr, 74.72.Dn 
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The charge stripe order [ljg] observed in the high- 
temperature cuprates superconductors is one of the 
most interesting current topics on the strongly corre- 
lated electron systems. In particular since the discov- 
ery of the coexistence with the superconductivity in 
Lai.6-a;Ndo.4Sr x Cu04 0, the mechanism of the stripe 
formation has been studied in many works. The numeri- 
cal study Q based on the density matrix renormalization 
group suggested that such a stripe phase can appear in 
the two-dimensional t-J model. On the other hand, the 
numerical diagonalization of the 4x4 t-J cluster with 
two holes || indicated that the stripe order occurs only 
in some low-lying excited states, rather than the ground 
state. The realization of the stripe order in the simple 
t-J model is still an open problem. 

It is well known that the t-J model should exhibit 
the phase separation for sufficiently large J/t. || The 
high temperature expansion suggested such a state is re- 
alized for J/t > 1. § Some small cluster calculations 
have shown that a larger cluster of the holes is stable 
rather than a pair even in more realistic parameter re- 
gion (J/t > 0.5). || In the present paper, we propose 
a possible mechanism of the stripe order formation due 
to the additional next-nearest-neighbor exchange interac- 
tion J' based on a naive argument valid in the phase sep- 
aration region of the t-J model. Since the next-nearest- 
neighbor hopping t' has been revealed to be quite large 
for Sr 2 Cu02Cl2 (f ~ 0.3t) J' is also expected to be 
finite in some real cuprates. Thus we consider the square- 
lattice t-t'-J-J' model to discuss on the mechanism of the 
stripe. We also calculate the three- and four-hole corre- 
lation functions of the 4x4 cluster with four holes, to 
examine the realization of the mechanism. 

We consider the two-dimensional t-J model in the pres- 
ence of the next-nearest-neighbor hopping t' and the ex- 
change interaction J'. The Hamiltonian is given by the 
form 
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H = -*E<i.j>, CT (4<T c ^ 

- f E<ij>', ff ( c L ci ^ + c L c s,°) 

+ J E<i,j>( s i - Sj - |wij) 



where E<ij> an d E<ij>' m ean the summation over 
all the nearest-neighbor and the next-nearest-neighbor 
sites, respectively. Throughout the paper, all the ener- 
gies are measured in units of t. We assume the next- 
nearest-neighbor exchange interaction is antiferromag- 
netic (J' > 0), as was revealed for La2Cu04 by the theo- 
retical study based on the ab initio calculation. JTo[ ] The 
antifcrromagnetic J' term can also be derived from the 
strong correlation expansion of the Hubbard Hamiltonian 
up to the order of t 4 /U 3 . jll] Since t' plays no essential 
roles in the following argument, we set t' = at first. 

Consider the naive argument to explain the hole pair- 
ing due to the antiferromagnetic short range order: a 
pair of holes sitting on the adjacent cites is more stable 
than two separated holes, because the former breaks 7 J 
bonds, while the latter 8 J bonds. Following the argu- 
ment, larger hole clusters are expected to be formed for 
sufficiently large J. In such a situation we consider the 
effect of J' . (We assume J' is not so large that the an- 
tiferromagnetic short range order is completely broken.) 
At first we compare the stability of three-hole cluster in 
two different shapes, shown in Figs. 1(a) and (b), respec- 
tively. The number of J bonds are the same between 
them, but (a) has one more broken J 1 bond than (b). 
When the antiferromagnetic short range correlation is de- 
veloped, the J bond should lead to the advantage of the 
energy, while the J' to the disadvantage, as far as J and 
J' are antiferromagnetic. Then (a) is expected to be more 
stable than (b). Thus the three hole cluster should prefer 
the line shape like (a) to the corner shape like (b). Next 
we consider the four- hole cluster with the two shapes, 
shown in Figs.l (c) and (d), respectively. In this case the 
number of J bonds is also different. One more J bond 
and two more J' bonds are broken in the shape (c) than 
(d). Assuming that the antiferromagnetic short range 
order is so large that the next-nearest-neighbor spin cor- 
relation is almost the same as the next one in amplitude, 
the line shape (c) is more preferable than (d) under the 
condition J' > J/2. This condition is easily revealed to 
be approximately valid in comparison between the line- 
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shaped and the square-shaped larger clusters with the 
same number of holes. Thus large line-shaped clusters 
of holes should be formed for sufficiently large J'. This 
naive argument is expected to give a possible mechanism 
of the charge stripe order. 



in the ground state of the finite cluster t-t'-J-J' model 
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FIG. 1. Schematic figures to discuss on the stability of the 
three-hole and four-hole clusters. 
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3-hole correlation functions 
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4-hole correlation functions 
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FIG. 2. Configurations of the many-hole correlation func- 



tions; (a)Cg\ (b)C^ s \ (c)C^ j and (d)C, 

In order to examine the realization of the mechanism 
of the charge stripe order discussed in the previous sec- 
tion, we calculate the three- and four-hole correlation 
functions defined as 



(3) 



-.(4) 
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(t' = 0). Cg^ and C^ t ' are supposed to represent a rela- 
tive strength of the stripe order, while Cpg and CpL mea- 
sure a tendency towards the ordinary phase separation. 
They are calculated for the 4x4 cluster with four holes, 
for which the ground state has the d-wave like rotational 
symmetry for J > 0.3. || (We neglect the other ground 
states which appear in smaller J regions for simplicity.) 
The calculated three- and four-hole correlation functions 
are plotted versus J' with fixed J (=0.6 and 0.8), in Figs. 
3 and 4, respectively. We detected a first-order transition 
(a level cross) at some critical value J' c (J' c depends on 
J.) and found that the line-shaped correlation is larger 
than the square-shaped one for J' > J' c , while it is re- 
versed for J' < J' c in both Figs. 3 and 4. It implies that 
the charge stripe order is possibly realized in the bulk 
system for sufficiently large J', in agreement with the 
mechanism proposed in the previous section. Then J' c is 
expected to be the boundary between the phase separa- 
tion and the stripe ordered phases in the thermodynamic 
limit. Plotting the calculated J' c for various values of J, 
we give a phase diagram in the J'- J plane for if = 
(solid circles) in Fig. 5. We can also understand that 
the excited state with the stripe order, which was found 
in the previous numerical study is stabilized by the 
next-nearest-neighbor exchange interaction in the upper 
phase in Fig. 5. 
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FIG. 3. Three-hole correlation functions versus J' with 
fixed J (=0.6 and 0.8). 
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FIG. 4. Four-hole correlation functions versus J' with 
fixed J(=0.6 and 0.8). 
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FIG. 5. Phase diagrams in J' -J plane for t'=0, -0.1 and 



0.1. 



The phase diagram for t' = in Fig. 5 indicates an 
interesting point that the stripe order is possibly real- 
ized even if J' is much smaller than J/2 in small- J re- 
gion around J ~ 0.4, which is realistic for the high-T c 
cuprates. Some recent theoretical analyses |l2|-|l4| on 
the simple t-J model actually revealed that the phase 
separation occurs even in such a realistic parameter re- 
gion. The present result of the phase separation-stripe 
boundary J c ~ 0.3 for J' = t' = in Fig. 5 is con- 
sisent with these results. It implies that the scenario of 
the stripe formation based on the next-nearest neighbor 
exchange interaction is possibly valid for real cuprates, 
although the presice phase boundary is still controver- 
sial. Note that the present anarysis does not distin- 
guish the static stripe order and the dynamical one like 
the charge strings, which was predicted by the phonon- 
induced polaron mechanism. |lq| It would be an inter- 
esting futre work to study on such a dynamical stripe, 
which may give some hints to explain the coexistence of 



the stripe order and the superconductivity observed in 
Lai. 6 -xNdo.4Sr a; Cu04. 

Finally, we consider the effect of the next-nearest- 
neighbor hopping t' in the present mechanism of the 
stripe formation. For this purpose, the phase boundaries 
between the stripe and the pairing (or phase separation) 
phases for £'=-0.1 (diamonds) and i'=0.1 (crosses) are 
shown in Fig. 5. The negative and positive t' are corre- 
sponding to hole and electron doping cases, respectively. 
The phase diagram suggests that the negative t' sup- 
presses the stripes, while the positive t' enhances it. The 
result agrees with the numerical studies fl6|,[l7j at least 
for small t', although they didn't consider J'. It implies 
that the stripe due to J' in the present mechanism has 
the same feature as the one which was investigated in 
those previous works. Actually Fig. 5 indicates that the 
stripe can occur even for J' = at least in the case of 
the positive t' . It would be more interesting to perform 
the same calculation for more realistic hole density, close 
to 1/8, if possible. (For example, the 32-site cluster with 
4holes is desirable, but it is difficult for the present com- 
puter systems.) 

The recent high-resolution inelastic neutron scatter- 
ing experiment JlT| indicated that the ring (four-spin) 
exchange interaction is more important to explain the 
observed spin- wave dispersion of La2Cu04, rather than 
the next-nearest-neighbor exchange interation. Thus we 
should also take the ring exchange interaction into ac- 
count for more quantitative study. 

In summary, we proposed a possible mechanism of 
the charge stripe formation based on the next-nearest- 
neighbor exchange interaction J' in the high T c cuprates. 
The many-hole correlation functions of the 4x4 lattice 
t-t'-J-J' model indicated that even small J' possibly in- 
duces the stripe order for realistic values of J. In addi- 
tion the next-nearest-neighbor hopping t' was revealed to 
suppress the stripe for t' < 0, but enhance it for t' > 0. 
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